\ 



L 

i 




• X-643-64-50 

TM X-55013 


SOLAR CYCLE CHANGES 

in 

INNER ZONE PROTONS 

N 64 27251 

| fa.t, &-£> 


MARCH 1964 



% 





OTS PRICE 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 


XEROX 

MICROFILM 


$ 5. 66 pk 

$ : 

f 


\ ' 


ABSTRACT ° 

Time dependent calculations of the inner Van Allen belt proton 
population show that large changes of the population up to a factor of 
50 will take place during the solar cycle. The effect is most pronounced 
for the region of B-L space corresponding to minimum altitudes of 300-700 
km. Because different energy protons respond to the changing solar cycle 
at different rates the proton energy spectra will change with time also. 





Solar Cycle Changes in Inner Zone Protons 


Robert C. Blanchard^ 
and 

Wilmot N. Hess 

The purpose of this paper is to calculate the expected changes in 
inner zone proton populations with time in the solar cycle. 

Freden and White [19591 identified the penetrating component of the 
inner zone of the Van Allen belt as due to energetic protons and measured 
the energy spectrum of the protons of E > 75 Mev. This and subsequent 
experimental work on Atlas rockets [Freden and White, i960; Freden and 
White, 1962; Armstrong, Harrison, Heckman and Rosen, 196l] has confirmed 
and extended this finding and we now have a well-established proton energy 
spectrum at L ^ 1,4, B ^ . 20 , which is about apogee for these flights, for 
near solar maximum (see Fig. l). The analysis of this spectrum [Freden and 
White, I960] has shown conclusively that the protons are produced by neutron 
decay. 

One of the features of this component of the radiation belt is the time 
constancy. Measurements after a solar flare [Freden and White, 1962] gave 
essentially the same fluxes of protons as before the flare. For L > 1.6 
changes in proton spectra are seen [Naugle and Kniffen, I96l] associated 
with solar proton events [Pizzella, 1961 ]. These may be explainable in 
terms of neutrons produced in the polar regions by solar protons [Lenchek, 
1962; Lingenfel ter , to be published]. For L < 1.6 only gradual changes in 
the proton populations have been observed on Explorer VII [Pizzella, 1 96 1 ] . 
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It was originally suggested [Pizzella, 1961 1 that such changes implied that 
neutrons could not produce all inner zone protons but this does not seem 
to be the case. Hess [ 1962] showed that slow changes in proton population 
were expected due to changes in the galactic cosmic ray flux during the 
solar cycle and more importantly due to changes in the upper atmospheric 
density during the solar cycle. We will now make a quantitative discussion 
of the expected changes in proton populations for L < 1.6. 

The continuity equation for protons is usually written [Singer, 1958; 
Hess, 1959; Freden and White, I960] 

= + S(E) - L(E) + ^ [ N(E) ~ v ] - 0 (1) 

and the steady state population of protons N(E) solved for. We are now 
interested in time variations so we will write the time dependent form of 
this equation 

= + S(E,t) - L(E,t) +“ [ N(E , t) ~ (t) v ] (2) 

The source term S(E,t) will use the neutron decay source strength of Hess 
for solar maximum S(E) multiplied by a function of time f(t) to consider the 
solar cycle variation in galactic cosmic ray flux as determined by McDonald 
and Webber [ 1961 3 
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S ( E t) = S(E) f(t) , ° F(t) (3) 

where f(t) varies from 1.0 at solar max to 1.25 at solar min in a similar 
way to Fig. 2. 

For the loss term L(E,t) following Freden and White [i960] we consider 
nuclear collisions 

L(E , t) = N(E,t) V £ a ± -Pi (t) (4) 

where *p^(t) is the average atmospheric density of atmospheric component 
i (0 2 , 0, N 2 or He) and <j± is the inelastic cross section for collisions 
(assumed geometric) . 

The last term in equation 2, the change in proton population due to 
slowing down, varies with the solar cycle due to the atmospheric density 
change 


dE r p (t) 1 dE- | 

dx 3 _2.69 x 10 _ dx 'OTP 


(5) 


where p (t) is the average atmospheric density of equivalent oxygen atoms. 
The variation of p £ (t) and p (t) with time in the solar cycle is the major 
reason for the change in proton population. 
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Average Atmospheric Densities 

A major part of this problem was to determine the average atmos- 
pheric densities used in calculating the loss rate of trapped protons by 
coulomb scattering. To do this the theoretical model of the upper atmos- 
phere developed by Harris and Priester [1962] was used. This model gives 
the time dependences of the atmosphere for both solar cycle and diurnal 
variations. It agrees well with several measurements of density by 
satellite drag (King-Hele, 1963; Bryant, to be published) and recently 

has been checked by preliminary data of densities from the Explorer XVII 
Spencer et al T 1 963 1 

satellite./ This model atmosphere is the most complete description of the 
time dependence of atmospheric densities available and agrees well with 
the current experimental data. We must perform several operations on the 
H and P model data to get it in form to use for this problem. Harris and 
Priester give the atom densities n^ of the several atmospheric constituents 
i in the form of , 


n i (h, t, S) 

where h is altitude above the earth, t is local time, and S is a model 
parameter related to, but not the same as, the intensity of the average 
10 cm solar flux, F, in watts/m 2 -cps x 10 22 . Recent studies of atmos- 
pheric densities [Harris and Priester, 1963] show that the model number 
S is the same as F near solar max but near solar min S < F. For example, 
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F = 70 corresponds to the model S =» 100. Fig. 2 shows the solar cycle 
variation of S we have used based on the last solar cycle. 

The first operation on the H and P atmosphere is to average over local 
time. The protons we will consider live long enough so that they will 
encounter the midday density bulge and the nighttime minimum many times and 
will average them out. 

Secondly, we perform a longitude average. As the particles drift 
around the earth their mirror points do not stay at constant altitude but 
rather follow a certain path in B-L space. Several such constant B-L rings 
are shown in altitude- longitude coordinates in Fig. 3* The particles dip 
lowest in the South Atlantic due to the nature of the earth's magnetic field. 
A B-L map of minimum altitudes is shown in Fig. 4- Due to the variation in 
altitude of the protons' mirror point, longitude averages were made using the 
B-L rings every ten degrees in both northern and southern hemispheres. This 
step gives the average mirror point density for a proton's motion. Walt 
has recently shown that the protons do not drift in longitude at a constant 
rate because of the variation of the magnetic field gradient and field line 
curvature. This effect is ignored here. It probably amounts to about a 
20 percent correction on average densities. Values of atmospheric densities 
are not known that well now anyway. 

The third operation on the H and P atmosphere is to average over a 
bounce from one mirror point to the conjugate mirror point. Since the 
protons live much longer than a bounce period, and since no change in direc- 
tion of the proton is considered during slowing, it is permissible to average 
the density in this fashion. A dipole field is assumed in carrying out this 


step. 



Ray [I960] gave the bounce averaging process as 


fp(x)ds 

p , 

I ds 

where the element of length along the particle’s orbit is 


ds 


vdt _ d£ 
cos a cos a 


where d £ is the element of length along the field line. Using 


d£ 2 = dr 2 4 d\ 2 


and 


r » L cos 2 X 


with the mirror equation and the expression for the magnetic field varia 
tion along a field line 


combining and substituting in equation 6 gives 


Jpi(x) A (x) dx 

'i * — 

J A(\) dX 


--7-- 


where 


. r g£S-. X 4±J_sin!_A j. \ fll ) 

[ cos 6 X (1 + 3 sin 2 x m )® - cos 6 \ m ( 1 + 3 sin 2 \)2] 

The subscript m here corresponds to the particles' mirror point. This 
bounce averaged density has been integrated on a computer. 

The last step is to combine the five i constituents N 2 , 0 2 , 0, He and 
H to give the averaged number of equivalent oxygen atoms 

8 p = 14 ~p (N 2 ) + P (0) + 16 p (0 a ) + 2 p (He) + p (H) ( 12) 

Values of ^ at different times in the solar cycle for L m 1.25 are shown in 
Figure 5. 

We can use the data on the rate of energy loss dE/dx for oxygen at 
NTP [Aaron, Hoffman and Williams, 1949] in equation 5 with this averaged 
density p to calculate the loss rate. 
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RESULTS 

Equation (2) is integrated stepwise starting with N(E) = 0 at t = 0 
to build up to an oscillating proton population which is the * same from one 

solar cycle to the next. In Figure 6 is shown the number of solar cycles 
after which the population change per cycle is less than 1$ for L = 1.25. 
After achieving this final condition the proton energy spectrum varies 
during one solar cycle as shown in Figures 7-11. The dotted curves top and 
bottom in these figures are what the proton spectrum would be if steady 
state conditions were achieved at solar max and solar min. Steady state 
clearly is not achieved for high energy protons or for high altitudes. 

It is also interesting that this calculation predicts a spectral peak 
above 100 Mev for certain conditions (Figure 7, year 3) due to the time 
lag in different energy protons adjusting to solar cycle changes. 

Figures 12 and 13 show time histories for various conditions. Large 
variations from solar max to solar min populations are seen in these 
figures. The amplitude of these changes for L = 1.25 and 1.40 are shown 
in Figures 14 and 15. 

A comparison of past experimental results and the calculations of 
this paper is shown in Table I. In discussing these experiments we can 
consider three types of information. 

In comparing time changes of the proton fluxes two situations occur: 

1. Some individual experiments run long enough (a year or more) 
so that they should see proton flux changes directly. From Table I, 



- 9 - 


experiments (c) and (d) are of this time. Although in experiment (d) 
it is not one experiment on one satellite, it is many identical experiments 
on essentially identical satellites, and therefore this experiment falls 
in this category. 

2. Different experiments on different satellites can be compared 
to provide data on proton flux changes. Data of this type is given in 
experiments (a), (e) and (f) in Table 1. 

The third type of data is given by: 

3. Direct measurement of the proton energy spectrum in one 
experiment. Experiments (a), (b) , (e) , (f) and (g) are of this type. 

In comparing flux time change experiments of type (l) with the 
calculations, the agreement is good. Two experiments covering the period 
1959 to 1962 give results consistent with the calculations. 

Comparing with type (2) data does not work as well. The observed 
changes of flux with time agree qualitatively than those predicted but 
not quantitatively. The experimental flux ratios are smaller than the 
calculated ones. Two remarks are in order here. First, type (2) data 
usually has larger errors attached to it than type (l) because it involves 

the systematic errors of two experiments while type (l) data involves no 
systematic errors--only statistical errors are involved in the time 
variations in the flux. Because of this the experimental flux ratios 
using type (2) data are probably not as accurate as those using type (l) 
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data. Secondly all the comparisons of type (2) data involve Explorer IV 
data; so there is no independent check of experimental consistency. 

We may have over-estimated the amplitude of the change of proton 
populations in this present calculation by as much as a factor of two 
because our solar cycle model uses S m 70 at solar minimum while recent 
data of Harris and Priester (1963) shows that S m 100 is probably more 
appropriate. Changing the shape of the model solar cycle (Figure 2) will 
also affect the results of the calculation. The solar cycle must rise 
more rapidly towards solar max. than it falls to solar min. in order to 
produce the changes in spectrum calculated here. Solar cycle 16 was 
quite symmetrical and should not give the calculated results. 

Of the five experiments of type (3) that help understand the proton 
energy spectrum the first two (a) and (b) in i960 showed a modest sized 
peak at about E 40 Mev. The more recent experiment (g) of Rowland et 
al in 1962 shown in Figure l6 does not show such a peak. The solid curve 
on Figure 16 compared with the data of Rowland is the results of the 
present calculations normalized. The agreement is fairly good. Recent 
experiments (f)(g) show a large low energy proton population but cannot 
be compared with early experiments because the early ones did not go down 
to such low energies. Figures 7-11 show that the calculated spectral shape 
should not change much during the period of these measurements- -from 6 to 
9 years solar cycle time. The proton fluxes should increase but the shape 
of the spectrum stays nearly the same. Striking changes in spectral shape 
should occur on the upswing of the solar cycle. The agreement with spectral 
shape in Figure 16 is not bad but the peaks in experiments (a) and (b) 
should not occur at t * 6. It may be that a neutron absorption as suggested 
by Freden and White (1962) is required to explain the peak. 
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All of the information presented so far on Figures 6 through 15 are 
mirror point fluxes M, that is, the flux of protons mirroring per unit 
volume at one particular B and L. For low altitudes where the atmosphere 
changes rapidly with altitude it is nearly correct to compare this flux 
with measured omnidirectional flux values* For these low altitudes the 
omnidirectional flux is very nearly the same as the mirror point flux. 

To show a more complete picture of the solar cycle proton changes we have 
converted to omnidirectional fluxes using equations (24) of Hess and 
Killeen (l96l). Figure 17 shows a calculated R-^ map of the 25 Mev proton 
omnidirectional fluxes J at solar max and solar min. An isoflux contour 
is clearly at lower R for solar minimum. 

In conclusion we have shown that large changes in proton fluxes will 
take place during the solar cycle for those regions of B-L space corres- 
ponding roughly to minimum altitudes of about 300-700 km. Changes in 
proton spectra will occur also. Comparison of the calculations with avail- 
able experimental information are not conclusive. Some kinds of experimental 
data agree quantitatively with the calculations. Other data agrees qualita- 
tively but not quantitatively. The crucial tests of the calculations will 
involve measuring changes in the proton flux and spectrum as we approach 
solar maximum during 1966-68. We know experimentally (King-Hele 1963) 
that the atmospheric density changes used in the calculation are reasonable. 
These calculated effects must take place with about the magnitudes , shown here 
unless there are features of the inner belt protons which we do not now 


understand. 
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FIGURE CAPTIONS 


1 . The experimentally measured inner zone proton energy spectrum (Freden and 
White 1962) for L ~ 1.4 B ~ .20. 

2. Solar cycle used in these calculations S is the Harris and Priester (1962) 
model parameter which is related to but not the same as the 10 cm flux 
intensity . 

3. Variation of the altitude of B-L rings with longitude. 

4. B-L map of mini-mum altitudes. 

5. Average atmospheric densities at different times of the solar cycle at 
different B-L points. 

6 . The number of solar cycles needed for a group of protons to come within 1 $ 

of its final oscillating population for L 3 1.25 starting with N(E) a 0 at 

t a 0 , 

7. Proton Energy Spectra at different times in the solar cycle for L » 1 . 40 > 

B s .225. The curves are labeled by the time in years from solar min. 

8 . Proton Energy Spectra at different times in the solar cycle for L » 1.40, 

B a .201. The curves are labeled by the time in years from solar min. 

9. Proton Energy Spectra at different times in the solar cycle for L a 1.40, 

B « .181. The curves are labeled by the time in years from solar min. 

10. Proton Energy Spectra at different times in the solar cycle for L s * 1.25> 

B 5 * .207. The curves are labeled by the time in years from solar min. 

11. Proton Energy Spectra at different times in the solar cycle for L =» 1.25> 

B =* .183. The curves are labeled by the time in years from solar min. 

12. Proton time histories for L * 1.40, B =* .225. 

13. Proton time histories for L se 1.40, B » .201. 

14 . The amplitude of the proton population change during the solar cycle for 

L = 1.25. 

15. The amplitude of the proton population change during the solar cycle for 

L = 1.4 and 1.6. 

16 The energy spectrum of inner zone protons measured by Rowland et. al . 1963 
* on an Atlas pod at L 3 1.27, B - .216 shown for comparison is the spectrum 
calculated in this paper for L m 1.27> B « .207, t a 9. 

17. An R-4 map of the omnidirectional flux of E » 25 Mev protons at solar max 

and solar min. 
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